Symmetries and weak (anti) localization of Dirac fermions in HgTe quantum wells 
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We perform a symmetry analysis of a 2D electron system in HgTe/HgCdTe quantum wells in the 
situation when the chemical potential is outside of the gap, so that the bulk of the quantum well 
is conducting. In order to investigate quantum transport properties of the system, we explore sym- 
metries of the low-energy Hamiltonian which is expressed in terms of two flavors of Dirac fermions, 
and physically important symmetry-breaking mechanisms. This allows us to predict emerging pat- 
terns of symmetry breaking that control the weak localization and antilocalization showing up in 
transverse- field magnetoresistance. 
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I. INTRODUCTION 

Solid state systems with massless Dirac charge car- 
riers have recently attracted an outstanding degree of 
attention of theoretical and experimental groups world- 
wide. This is mainly due to the impressive experimental 
progress in two physically related directions of research: 
graphene 1 ^ and topological insulators^. 

The notion of a topological insulator (TI) refers to a 
bulk insulator with gapless surface states occurring due 
to topological reasons. The simplest example of a topo- 
logical insulator is the 2D electron gas in a strong mag- 
netic field. At the quantum Hall plateau, the gap between 
Landau levels in the bulk is penetrated by a fixed inte- 
ger number of chiral edge states providing the quantized 
value of the Hall conductance. The integer quantum Hall 
edge is thus a topologically protected ID conductor real- 
izing the group Z. 

A novel class of TIs2r— requires strong spin-orbit inter- 
action in the absence of magnetic field (i.e. it is realized 
in systems with preserved time-reversal invariance). This 
type of TIs was discovered in HgTe/HgCdTe structures 
by the Wurzburg group£. Strong spin-orbit interaction in 
HgTe leads to the inverted band gap in this semiconduc- 
tor. As a result, the electron and hole bands are cross- 
ing near the boundary of the sample giving rise to the 
two counter-propagating helical edge modes. The time- 
reversal symmetry of the system leads to the topological 
protection of these edge modes. 

Voltage applied to such a sample results in the appear- 
ance of the perpendicular spin current. This phenomenon 
is known as the quantum spin-Hall effect (QSHE). The 
robustness of the effect with respect to disorder makes 
it an extremely promising tool for applications. As 
a simplest example, the conversion between the usual 
charge current and spin current occurring in QSHE 
can be used for generation and detection of spin cur- 
rents. The existence of a non-localized conducting chan- 
nel at the edge of an (appropriately manufactured) two- 
dimensional HgTe/HgCdTe quantum well was experi- 



mentally demonstrated in Refs. I^flfjl These experiments, 
showing that HgTe/HgCdTe quantum wells provide a re- 
alization of a novel remarkable class of materials — Z2 
topological insulators — opened a new exciting research 
direction. 

A 3D realization of a Z2 TI was discovered in the ex- 
periment by Princeton group^ where crystals of Bii-^Shr 
were studied (later other Bismuth compounds such as 
BiTe and BiSe were also shown to be 3D topological insu- 
lators). The 3D topological insulators exhibit very strong 
spin-orbit interaction also leading to the band gap inver- 
sion. This results in the appearance of gapless states on 
the surface of the sample forming a 2D topologically pro- 
tected metal. The dynamics of the surface states is gov- 
erned by the same massless Dirac Hamiltonian that has 
previously appeared in graphene. The main difference 
between graphene and the surface of a 3D topological in- 
sulator is the lack of the spin and valley degeneracy in 
the latter case. A combined effect of topology, disorder 
and interaction on a surface of 3D topological insulators 
was demonstrated to establish a novel critical 2D state 
with the low-temperature conductivity of order unityJ^ 

A similar critical state is expected^ to separate the 
normal and topological insulator states in a 2D Zi topo- 
logical insulators. In particular, such critical state should 
occur exactly at the QSHE transition in HgTe/HgCdTe 
quantum wells of critical thickness, when the 2D bulk 
gap is tuned to zero. The properties of this topologically 
protected metal reflect the quasirelativistic Dirac nature 
of carriers, similarly to the interference phenomena in 
graphene 

A two-dimensional metallic state in HgTe/HgCdTe 
quantum wells can also be realized in the presence of 
the bulk gap (either on the TI on the normal insula- 
tor side) by shifting the chemical potential (with a help 
of the gate) away from the gap. Such two-dimensional 
TI away from the TI regime represents a 2D spin-orbit 
metal, whose properties (in particular, interference cor- 
rections to the conductivity and low-field magnetoresis- 
tance) may still reflect the Dirac-fermion nature of car- 
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riers. In recent transport experiments the magnetoresis- 
tance of bulk-conducting HgTe/HgCdTe quantum wells 
was studied, both for inverted (thick wells )iSr— and nor- 
mal (thin wells) 20 i 21 band structures. 

Although some aspects of the interference corrections 
in TI systems away from the TI regime have been al- 
ready addressed theoretically-^ there is a clear need in 
a systematic theory describing the whole variety of ex- 
perimentally accessible regimes. The purpose of this 
paper is to develop such a theory. Analyzing the un- 
derlying Dirac-type Hamiltonian and physically relevant 
symmetry-breaking terms, we identify parameter regimes 
with different symmetries. The corresponding symmetry- 
breaking patterns determine the form of quantum in- 
terference magnetoresistance, i.e. the weak localization 
(WL) vs weak antilocalization (WAL) and corresponding 
prefactors. 

The paper is organized as follows. In Sec. [TT] we ana- 
lyze symmetry properties of the Bernewig-Hughes-Zhang 
Hamiltonian and various physically relevant symmetry- 
breaking perturbations. This allows us to establish 
emerging patterns of symmetry breaking. In Sec. Mil we 
use these results to calculate the interference corrections 
(weak localization and weak antilocalization magnetore- 
sistance) in various parts of the parameter space. Section 
IIVI contains a summary of our results and a brief com- 
parison to available experimental data. 



II. EFFECTIVE HAMILTONIAN AND 
SYMMETRIES 

The standard approach to the description of electronic 
states in semiconductors is based on the k • p method. 
This method assumes that exact wave functions in the 
periodic potential of the crystal are known at some given 
point in the Brillouin zone. Normally, the states with 
zero quasimomentum (T point) are used. Then the ef- 
fective Hamiltonian is constructed in the basis of these 
states by expansion in small momentum k. The k ■ p per- 
turbation theory thus provides the model Hamiltonian 
for the states close to the T point in terms of few matrix 
elements calculated with respect to exact eigenstates at 
p = 0. The symmetries of the underlying system impose 
a number of constraints on the structure of the effective 
theory. In fact, a general form of the effective 2D Hamil- 
tonian can be derived from the symmetry grounds, see 
Ref. [T3I without invoking the microscopic (tight-binding 
or Kane-model) Hamiltonian. 

The only intrinsic symmetry of the 3D Hamiltonian is 
time-reversal (TR): H = s y H*s y . In the clean case, the 
point symmetry group is Td (tetrahedral) . Spatial inver- 
sion interchanges Hg and Te sublattices of the crystal. 
Neglecting the difference between Hg and Te, the point 
group becomes Oh (cubic). The approximate inversion 
symmetry leads to the block-diagonal structure of the 
BHZ Hamiltonian. 



A. BHZ Hamiltonian 

Effective Hamiltonian for a narrow symmetric HgTe 
quantum well (QW) was derived in Ref. [5| in the frame- 
work of the k • p method. The Bernevig-Hughes-Zhang 
(BHZ) Hamiltonian has a 4 x 4 matrix structure in the 
spin (sign of the ^-projection of the total momentum J, 
where the z-axis is perpendicular to the QW plane) and 
El - HI space, 

h - ( h( ® ^ m 
HV = (fl + m t] A £ x + l %\)- ( 2 ) 

v ; \A{k x - iky) e(k)-m(k)y v ; 

Here we have used the form given in Refs. I4l fl2lfl3l with 
the following arrangement of components in the spinor: 
E1+, H 1+, El — , HI — . Introducing the Pauli matrices 
<7a,x,y,z f° r the El — HI space and so.x.y.z for the up — 
down (spin) space (with ctq and sq unity matrices), the 
effective Hamiltonian can be written as 

Hbhz = e(k)a s +mQs.)a z s + Ak x a x s z - Ak y a y s . (3) 

The functions e(k) and m(k) are effective energy and 
mass. Within the k • p expansion in the vicinity of the T 
point, they are 

e(k) = C + L>k 2 , m(k) = M + Bk 2 . (4) 

The two phases of normal and topological insulator corre- 
spond to M > and M < 0, respectively. The sign of M 
changes at the critical thickness d c of the QW of about 6.2 
nmi Parameters A, B, and D are positive with B > D. 
The parameter C includes chemical potential and can be 
varied by changing the electron concentration with the 
gate voltage. 

The Hamiltonian i?BHZ breaks up into two blocks act- 
ing independently in the spin-up and spin-down sub- 
spaces with the spectrum 

E ± (k) = e(k) ± ^A 2 k 2 +m 2 (k). (5) 

The standard wayi^ to introduce disorder in the model 
is to add the fully diagonal term 

H dis = V(r)a s (6) 

with a random potential V(r) to the effective Hamilto- 
nian Hbhz- This corresponds to a relatively smooth (on 
the scale of the quantum well thickness) disorder. In 
particular, this model correctly describes charged impu- 
rities located at a certain (sufficiently large) distance R 
from the QW, e.g., in the doping layer. In this situation, 
within the QW the impurity potential is almost constant 
in z direction (across the QW), so that it does not break 
the z — > — z symmetry of the QW. Alternatively, one can 
consider, e.g., placing short-range impurities exactly in 
the middle of the QW: such impurities would also pre- 
serve the z — > — z symmetry and hence would not give 
rise to the mixing of the up and down blocks in Hbhz- 
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In general, the Hamiltonian Hbhz + #dis possesses two 
symmetries. Apart from the physical time-reversal sym- 
metry .ff(k) = SyH*(—\i)Sy, that relates the two spin 
blocks (Kramers doublets), the Hamiltonian commutes 
with s z . This allows one to consider the two blocks sep- 
arately. An extra symmetry can arise at some specific 
values of energy. In particular, h(h) acquires its own 
time-reversal symmetry when the mass m(fcj?) is zero. 
This happens in the inverted regime (M < 0) when 
M + Bkp = 0. 

It is also possible to achieve approximate orthogo- 
nal symmetry close to the bottom of the band, when 
| Af | > {Ak F , Bkp}, and at very high energies Bk 2 F » 
{|M|, Akp}. In the absence of block mixing, this would 
lead to WL. On the other hand, the regions of WL be- 
havior at stronger magnetic fields (crossing over to WAL 
in weaker fields) were found in an experiment on simi- 
lar structure d 19 ' 20 , indicating that one can indeed achieve 
the regime of (approximate) orthogonal symmetry in this 
class of devices. 

Thus we have identified three possible relevant sym- 
metries of the model: physical time-reversal, spin, and 
block-wise time-reversal (either "symplectic" or "orthog- 
onal" ) . The latter symmetry can occur only in the pres- 
ence of the spin symmetry, when the Hamiltonian breaks 
into two blocks. When all three symmetries are present, 
the two copies of symplectic or orthogonal class are real- 
ized leading to double WAL (i.e., with magnitude twice 
larger than the usual WAL correction) or double WL cor- 
rection (with the same prefactor as in the conventional 
spinful orthogonal class). If the block- wise time-reversal 
symmetry is broken while spin symmetry is preserved, 
we obtain two copies of the unitary class with no inter- 
ference corrections (in the lowest one- loop order) . If the 
spin symmetry is also broken, the system realizes a single 
copy of symplectic class yielding the usual (not double) 
WAL correction. Finally, if the physical time-reversal 
symmetry is broken (e.g., by magnetic impurities), the 
single copy of unitary class is realized with no interfer- 
ence corrections (in the lowest order). 



B. Symmetry breaking mechanisms in 2D 

Let us list possible mechanisms of symmetry breaking 
in a 2D quantum well of HgTe. 

Block-wise time-reversal symmetry is not exact from 
the very beginning. It is exact for massless Dirac 
fermions but is violated by the presence of m(k) =/= 0. 
In Hbhz, this symmetry occurs only at one particular 
electron concentration. Since disorder mixes the states 
in the energy window of order 1 Jt near Fermi surface, it 
will inevitably break this symmetry. Another, and more 
relevant, possibility is the detuning of the average elec- 
tron concentration from the point M + Bkp = 0. 

Spin symmetry is broken by one of the following mech- 
anisms: 

(i) The block-diagonal structure of Hbhz is not exact. 



Off-diagonal elements in the effective Hamiltonian arise 
due to bulk inversion asymmetry (BIA) of the zinc-blende 
lattice of HgTe and have the formic 
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+ S e (a Q + a z )(k x s x - k y s y ) + 5 h (a - o z ){k x s x + k y s y ). 

(7) 

Here the first term with Ao comes from the fc 2 matrix el- 
ements connecting electrons from the Tq band and heavy 
holes from the Tg band that have opposite spins (such 
coupling is absent in a spherically symmetric Kane model 
in a bulk system) J£ This term leads to the splitting of 
the spectrum into four branches: 



E ± (k) = e(k) ± v/(A|fc|±A ) 2 +m 2 (k). 



(8) 



The linear-in-fc terms with S e and Sh arise after the pro- 
jection of the bulk cubic Dresselhaus terms for electrons 
and holes, respectively, onto the QW. Typically, the spin- 
orbit interaction for holes is stronger than for electrons: 
Sh fie- The role of Dresselhaus- type terms increases 
with increasing kp. These BIA mechanisms of break- 
ing s z symmetry are characterized by the corresponding 
symmetry breaking rates, 

1 



1 

e.h 



A6V, 



{5 e ,hkF) 2 T, 



(9) 
(10) 



within the Dyakonov-Perel spin relaxation mechanism 
(here r is the transport scattering time due to disorder). 

The BIA-induced terms are often neglected in the con- 
sideration of the quantum spin Hall effect in view of small 
coupling constants. Nevertheless, the break-down of the 
spin symmetry certainly occurs at sufficiently long scales 
(low temperatures) , leading in the infrared limit to a sin- 
gle copy of a symplectic-class system (WAL without dou- 
bling). 

(ii) Another way to break s z symmetry is provided by 
the Rashba term arising in an asymmetric well (e.g., due 
to a finite gate voltage); 12 ' 13 



H R = 



'y + kySx). (11) 

Here we keep only the leading (linear-in-fc) Rashba 
term connecting El up and El down subbands; other 
terms involving HI bands are of higher order in k. 
Within the Dyakonov-Perel mechanism, the correspond- 
ing symmetry-breaking rate is 

— ~ (r k F ) 2 T, (12) 

TR 
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(iii) The diagonal structure of Hbrz is broken by atom- 
ically sharp impurities or interface roughness. This type 
of disorder, combined with the strong spin-orbit cou- 
pling in HgTe, leads to random non-diagonal terms in the 
Hamiltonian. Furthermore, the s z symmetry can also be 
broken by finite-range (smooth on the lattice scale) im- 
purities located in the quantum well. Such impurities 
break locally the z — > — z symmetry of the QW and can 
be thought of as a kind of random Rashba term. Near 
the transition point, the leading (linear-in-fc) contribu- 
tion comes from the disorder-induced mixing of El-up 
and El-down blocks, similarly to the Rashba term. When 
disorder is created by distant charged impurities, such 
mixing is suppressed by an additional factor d/R, where 
d is the QW thickness and R is the spacer width. 



C. Symmetry breaking mechanisms at the 
boundary 

When the dephasing length is longer than one of the 
dimensions of a sample, the boundary conditions become 
important for the interference effects. This is, in par- 
ticular, the case for universal conductance fluctuations 
in a coherent sample. Another example is the conduc- 
tance of narrow stripes of HgTe in the regime L^, 3> W. 
While the width of the quantum well is of the order of 
few nanometers, its lateral width W can be restricted to 
about hundred nanometers. In such a restricted quasi-lD 
geometry, boundary conditions play a crucial role. Let 
us analyze the symmetries of the BHZ Hamiltonian with 
the boundary conditions. 

We first assume that the boundary does not violate the 
spin symmetry and the Hamiltonian decomposes into two 
time reversed blocks h(k) and h*(— k). Then the block- 
wise time reversal symmetry of h(k) will be inevitably 
broken at the boundary. The easiest way to see this is 
to consider the Dirac limit of the Hamiltonian neglecting 
the mass in Eq. ([3]). Schrodinger equation with such a 
linear-in-fc Hamiltonian is a system of two coupled lin- 
ear first order differential equations. Boundary condi- 
tion for such a system should be of the zero order in 
momentum, i.e., some linear constraint on the compo- 
nents of the wave function is imposed at the boundary. 
In a general form, we can write b T ip = with some two 
component spinor b. The elements of b may depend on 
the conserved momentum parallel to the boundary. The 
time-reversal symmetry tells us that, if ijj is an eigenstate, 
then Syip* is another eigenstate with the same energy. If 
boundary conditions preserve this symmetry then the lat- 
ter eigenf unction would obey b T s y ip* — 0. Equivalently, 
bhyip = 0. Thus we have two linear conditions on tp 
instead of one. These conditions are consistent only if 
the vectors b and s v b* are linearly related: b = as y b* 
with some constant a. Multiplying the last expression by 
&t from the left, we obtain b^b — ab^Syb* — atr(s y b*tf). 
The left-hand side is strictly positive while the right-hand 
side contains a trace of the product of antisymmetric (s y ) 



and symmetric (b^b*) matrices and is hence zero. This 
apparent controversy proves that any boundary condi- 
tions for a single copy of massless Dirac Hamiltonian will 
inevitably break its time-reversal invariance. 

The simplest way to see this is as follows: in order to 
produce a "wall" for a single Dirac-fermion species, one 
has to open a gap by switching on a big mass near the 
boundary, which breaks the effective time reversal sym- 
metry already after a single boundary scattering event. 
This fact is well known in the context of graphene studies, 
where in the absence of the intervalley scattering it is im- 
possible to confine Dirac quasiparticles without opening 
the gap at the boundary. 

An alternative boundary condition rp = for the BHZ 
Hamiltonian is widely used in literature. In order to ap- 
ply this boundary condition, the quadratic terms in the 
Hamiltonian should be retained. The relation ip = is 
invariant under time-reversal transformation. Neverthe- 
less, the time-reversal symmetry is broken at the bound- 
ary even in this case. A detailed proof of the symmetry 
breaking is relegated to Appendix [AJ 

We thus conclude that the block-wise symplectic time- 
reversal symmetry is always broken near the boundary. 
In particular, in a quasi- ID system the corresponding 
symmetry-breaking rate is given by 



where W is the stripe width. This r e< ige is just the average 
time for an electron to diffuse across the system and get 
aware of the boundary conditions. 

When the block-wise time-reversal symmetry is of the 
orthogonal type (e.g., when the Fermi energy is located 
near the bottom of the spectrum), the boundary con- 
dition ip = does not introduce additional symmetry 
breaking, in contrast to the "Dirac case" . Indeed, the sys- 
tem belonging to the conventional orthogonal class can 
be confined by large potential. Then the "relativistic cor- 
rections" would interplay with the boundary scattering 
just in the same way as with the impurity scattering. 

So far, we have considered the boundary conditions 
preserving the spin symmetry. The spin symmetry may 
be broken by the edges of a 2D sample if the edges are not 
ideal in z-direction. This situation, which can be modeled 
by short-range impurities located near the boundaries, 
seems to be quite likely in a realistic setup. In this case 
the only remaining symmetry is the physical (symplectic) 
time-reversal symmetry. 

III. INTERFERENCE CORRECTIONS FROM 
SYMMETRY CONSIDERATION 

In this section we describe the general symmetry-based 
formalism for calculating the interference corrections in 
an infinite 2D sample. We closely follow the approach 16 
developed for Dirac fermions in a disordered graphene 
(see also Ref. [H). 
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Quite generally, the existence of a singular (logarith- 
mic in T or B) conductivity correction is related to the 
presence of a certain time-reversal (TR) symmetry that 
acts on an operator O according to 

T : Oh) U^ 1 O t U, (14) 

where U is some unitary operator (note that the momen- 
tum operator changes sign under transposition) . For the 
4x4 matrix Hamiltonian one can introduce 16 possible 
TR symmetries 

Tij : O h> a i s j O T s j a u (15) 

that would correspond to 16 soft modes (Cooperons). 
For example, the Hamiltonian H20 = i e (k) + V(r)}aoSQ 
(physically, this Hamiltonian corresponds to spinful elec- 
trons in two valleys of a normal metal, without any spin 
and valley mixing effects) is invariant under all these 16 
symmetries. This gives rise to 16 logarithmic corrections 
to the conductivity of the form 

2 

Saij = -c„ ^ In (J^j , dj = ±1. (16) 

The sign factors cy here are determined by the sign of the 
time-reversal operator squared T 2 = ±1. In the present 
case we have T?- = (<JiSj)(aiSj)* . In particular, 

a s (a s a )* = 1 =>■ T 2 = 1, (17) 
<j y s (a y s y = -1 =>■ T 2 = -1, (18) 
(T y s y (a y s y y = 1 =>> T^ y = 1, (19) 

and so on. One can see that out of 16 TR symmetries only 
6 (namely, T 0y , T xy , T y0 , T yx , T yz , T zy ) contribute with 
the minus sign, yielding for H20 the total conductivity 
correction 




This is the WL correction for two independent copies of 
orthogonal symmetry class. 

In general, not all the possible TR symmetries are 
respected by the dominant term in the Hamiltonian. 
Therefore, as a first step in calculating the logarithmic 
correction one should retain only the TR symmetries of 
the dominant term, setting all other = 0. The sub- 
leading terms in the Hamiltonian may further break the 
remaining TR symmetries on the scale r s b, introducing 
gaps in the soft modes and thus cutting off the logarith- 
mic terms by T 8 b- As a result, in the low-T limit, the 
singular term ln^/r) remains only for those TR sym- 
metries that are preserved by all terms in the Hamilto- 
nian. With increasing T, becomes shorter than the 
symmetry breaking time r s & and the system crosses over 
to another symmetry class. 



A. Weak antilocalization in HgTe quantum wells: 

E F > m(k F ) 

Let us now return to the case of HgTe and employ this 
machinery to it. In order to directly use the results of 
Ref. [IB, we interchange El and HI states in the spin- 
up block. In this basis (H1+,E1+,E1-,H1-), the BHZ 
Hamiltonian © takes the form: 

Hbhz = e(k)<r So + [-m(k)a z + Aka]s z . (21) 

The linear-in-fc (massless Dirac) term in this 4x4 Hamil- 
tonian has the same structure as in Ref. [IH (there the 
matrices Tj corresponding to two graphene valleys played 
a role of Si here). 

We will now consider the situation when this term in 
the Hamiltonian is dominant, i.e. when the Fermi energy 
Ep 3> m(fcir) is in the range of almost linear spectrum 
E± ~ ±A|fc|, so that the mass term as well as the spin- 
symmetry breaking terms can be considered as perturba- 
tions. In this representation, the s^-symmetry breaking 
terms read: 

a y )s x (22) 

+ 5-(k X <Jy - ky(J X )S yi (23) 

Hr = r [(k x a y + k y a x )s x - (k x a x - k y a y )s y ], (24) 

where S± — 5h± 5 e . 

The massless Dirac Hamiltonian ("graphene Hamilto- 
nian" ) 

H A = A(k x a x + k y a y )s z (25) 
is invariant under the following four TR symmetries: 



T • 


O 1 y a x s x O T a x s x , 


r ~r2 

XX 


= 1, 


(26) 


T y0 : 


O 1 ^ <7yS Q O T <7 y S , 


rji'2 

1 y0 


= -1, 


(27) 


T ■ 


O 1 ^ (T y s z O T a y s z , 


1 yz 


= -1, 


(28) 


T ■ 


O 1 ^ a x s y O T a x Sy, 


rpl 

1 xy 


= -1, 


(29) 



which in Ref. were denoted as To,T x ,T y ,T z , respec- 
tively. As a result, the conductivity correction for two 
independent copies of massless Dirac fermions is positive: 

fc = -(l-3)-^-ln(^0 =2x ^rki(^) , (30) 

which is a doubled WAL (two copies of symplectic class). 

Let us now take into account the mass term Hm = 
— m(k)cr 2 s 2 in Hbhz- This term is invariant under T xx 
and T xy while breaking down T v q and T yz . The corre- 
sponding symmetry breaking rate l/r m was calculated 
in Ref. 14. The conductivity correction in this case can 
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be written as 

.2 
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\ T <t> 



2x ^r ln - + — • ( 31 ) 



At lowest T, when t$ 3> r m , we have 



2 x —-In 

2irh 



(32) 



which is not singular in T. This corresponds to two copies 
of a unitary class. With increasing T the symmetry- 
breaking term l/r m becomes unimportant and the sys- 
tem crosses over to two copies of symplectic class, Eq. 

(3D). 

So far, we have considered the block-diagonal s z - 
symmetric Hamiltonian -Hbhz • Let us now take into 
account the inversion-asymmetry terms Hbia and Hr. 
The main difference between HgTe and graphene is the 
structure of such symmetry breaking terms. The inter- 
valley disorder scattering in graphene preserves only T xx 
symmetry. As a result, when the valley mixing is strong, 
the quantum correction to the conductivity is negative 
(WL): graphene with mixed valleys belongs to the or- 
thogonal symmetry class. In contrast, the mixing of the 
up/down blocks in HgTe QWs involves an additional spin 
structure which generically kills the orthogonal T xx sym- 
metry. Generically, the only remaining TR symmetry is 
the symplectic T xy symmetry surviving the SO interac- 
tion. Therefore, at lowest temperatures the HgTe QW is 
expected to be a single copy of a symplectic system with 
ordinary WAL. 

The k- independent term Ao<7 z s x in Hbia violates T xx 
and T y o, while preserving T yz and T xy . This leads to 
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(33) 



Two situations are possible: (i) ta 3> r m and (ii) ta <C 
T m . Starting at high T, with lowering T in case (i) the 
system first crosses over from two copies of symplectic 
class, Eq. (I3TJ1) . to two copies of unitary class, Eq. (f3"2"j). 
and then at ~ ta to a single copy of a symplectic class 
with ordinary WAL: 



So 



2wh \t J 



We denote such an evolution as 2Sp 
case (ii), for ta -C we get 



2U 



So 
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' m J 


T 4>\ 



(34) 
ISp. In 

(35) 



which corresponds to the crossover 2Sp — > ISp at ~ r m 
(2WAL -» 1WAL). Note that at t ~ t a the number 
of soft modes does not change, despite the transitions 
between the spin blocks. 

Finally, in the presence of Dresselhaus and / or Rashba 
terms, the only effective TR symmetry is T xy . The effect 
of short-range impurities within the QW is similar to 
that of the Rashba term. Denoting the total symmetry- 
breaking rate due to such terms as 1/tso, we get (here 
neglecting the A -term): 



5a = 



2-Kh 



mil- 

T 4> 



2 In 



r 

tso 
r 

h 

T 4> 



T 
Tm 



T 

TSO 



In 



(36) 



Again, for r m <C rso, we have 2Sp 2U — > ISp, whereas 
for T rn 3> rso the mass term is not important (2Sp — ¥ 
ISp). 

The general formula for the conductivity correction 
at Ep 3> m(kp) including all symmetry-breaking terms 
combines Eqs. (|3"3")l and (j3"rJj) : 



5a 



e 

2nh 



In 



T 



T 



T 
TSO 



, I T T T T \ 

In ( — + — + — + 

S<P T m T A TSO/ 



-In 



T 
T 4 , 



T 
7~ri i 



-L -hi 

TSO / T0 



(37) 



Here is the dephasing rate due to inelastic processes, 
r is elastic transport scattering time, ta was defined in 
Eq. © , 1 /tso is the total spin-orbit rate describing the 
combined effect of Dresselhaus and Rashba terms, Eqs. 
(fTT)|) and (TT2")) (as well as the SO impurity scattering), 
1 /r m is the rate of breaking effective time- reversal sym- 
metry within each spin block due to finite mass of Dirac 
fcrmions. This mechanism of symmetry breaking is anal- 
ogous to the one in the case of interplay between Rashba 
and Zeeman splitting: the Rashba term tends to fix the 
spin direction in the 2D plane according to the momen- 
tum of a particle, whereas the Zeeman field tends to align 
the spin in the direction of magnetic field (out of the 2D 
plane). Without impurity scattering, each state is char- 
acterized by a certain direction of spin. The symmetry- 
breaking rate can be estimated using an analogy with 
the Dyakonov-Perel mechanism, t^ 1 ~ A^t, where ef- 
fective spin-precession frequency A m is only non-zero in 
the presence of disorder: 



m(kp) 1 
Ak F t : 



which yields 



m(kp) 
E F 



(38) 



(39) 
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in the regime of interest. A rigorous calculation of a 
Cooperonil (which can be simplified by dressing the im- 
purity potential by Dirac factors, thus reducing the prob- 
lem to a spinless one in the presence of a peculiar disorder 
potential) confirms this estimate. This result agrees with 
Ref. in this limit. 

Assuming for simplicity an energy-independent disor- 
der scattering rate (which, in fact, is jrenerically not the 
case for Dirac particles, see, e.g., Ref. [la), we can express 
the symmetry-breaking length l m = (flTm) 1 / 2 in terms 
of the carrier density n: 



.4 



\/n 



(40) 



where the sign ± distinguishes between the inverted (— ) 
and normal (+) band structures. 

We have identified the only true soft mode (corre- 
sponding to the true Kramers TR symmetry) , which gov- 
erns the localization properties at lowest temperatures. 
It gives rise to a single WAL correction (single symplec- 
tic class: ISp). This is not surprising, as the symplectic 
TR symmetry is the only true TR symmetry in a system 
with strong spin-orbit interaction. On the other hand, 
on short time scales (equivalently, high temperatures) the 
symmetry may be higher if the corresponding symmetry- 
breaking terms are not yet operative. Depending on the 
values of symmetry-breaking rates, there may be differ- 
ent patterns of symmetry breaking: 

(i) 2Sp — > 2U — > ISp. This first scenario is realized 
when the mass term in Hbhz is more important 
than the s 2 -symmetry-breaking terms: r <C r m <C 
min[rA, T50]. In this case, at high temperatures 
the system behaves as two copies of symplectic class 
(2Sp) with doubled WAL. In the intermediate range 
of temperatures r m <C <C min[rA, tso], the sys- 
tem behaves as two copies of unitary class (2U) 
with no T-dependent interference correction. Fi- 
nally, at ^ min[rA, tso] the two spin blocks 
are completely mixed and the system reaches its 
generic spin-orbit state ISp (single copy of sym- 
plectic class with ordinary WAL). 

(ii) 2Sp — > ISp. This scenario is realized in the 
case of the following hierarchy of scales: r <C 
min[rA, tso] T m, when the block mixing is faster 
than the breaking of the effective TR symmetry 
within each block. As a result, there is no room for 
the unitary class in this case. 

(iii) ISp. When the complete mixing of spin blocks 
is very fast, which happens when rgo ^ T or 
max[r m ,TA] < r, the system behaves as a single 
copy of symplectic class (ISp) in the whole diffu- 
sive regime (ordinary WAL). 

Let us estimate which of these scenarios can be ex- 
pected in experiments on HgTe QWs. According to Ta- 
ble I of Ref. 0, the value of the BIA-induced splitting 



A in HgTe QW near the QSHE transition is A ~ 
0.0015 - 0.002eV - 15 - 20K. When the 2D mean free 
path is 0.1 — 0.5/im, one finds (using A ~ 3.65 — 3.9 eV 

A) 



A r ~ 0.5 - 2.5 ~ 1. 



(41) 



This implies that for such values of the mean free path 
there is no room for scenario (i). The conductivity cor- 
rection for Aqt > 1 has the form 



5a 



e 

2nh 



(42) 



In 1 h - In — 

\T<f, T m T S J T$_ 

Then for min[r m , tso] 3> T the pattern (ii) 2Sp — > ISp is 
realized, while for min[r m ,Tso] ^$ t there is single copy 
of symplectic system (ISp) in the diffusive regime. 

In a perpendicular magnetic field, a magnetoresistance 
arises due to the suppression of the interference correc- 
tion by magnetic field. In the diffusive regime L$, Lh 3> I 
[where Lh = (he/ eB) 1 / 2 is the magnetic length and I is 
the mean free path], the effect of magnetic field can be 
described by the replacement 

1 1 



D 

H 



(43) 

in the above zeio-B formulas (here D is the diffusion 
coefficient). A more accurate expression (valid also in 
the limit of weak magnetic field, Dt^,/L h -C 1) has the 
standard Hikami-Larkin-NagaokaF^ form with digamma 
functions replacing logarithms. For a single copy of 
symplectic system, Eq. (|34[) . the magnetoconductivity 
Aa(B) = a(B) - a{0) has the from 



e 



MB) = -ir-r 



HI-,— 

2nh \ T<j, ' B tr 



x 



I) 



-In a;, (44) 



where B tr = h/(2el 2 r ) is the value of magnetic field for 
which Lh is equal to the transport mean-free path l tr and 
ip(x) is a digamma function. When the zevo-B expression 
contains a combination of several logarithms represent- 
ing the contributions of various channels, in the formula 
for the magnetoconductivity, each of logarithms is to be 
replaced according to Eq. (l44l) . with 1/r^ — > 1/t^ + 1/t s („ 
where 1 /r s b is the total symmetry breaking rate for given 
channel. Weak antilocalization leads then to a positive 
magnetoresistance. 

Importantly, in the regime Ep 3> m(kp) considered 
above there is no room for weak localization behavior 
which would manifest itself as a negative magnetoresis- 
tance. However, in Refs. [I9II20I a crossover form positive 
to negative magnetoresistance was observed with increas- 
ing magnetic field in samples with not too high carrier 
densities. Such a crossover is characteristic to conven- 
tional spin-orbit systems with tso 3> T - In order to 
describe this behavior in HgTe quantum wells we have 
to consider the opposite limit Ep <C m(kp), where the 
off-diagonal ( "spin-orbit" ) terms in the Hamiltonian can 
be considered as a small perturbation. 
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B. Weak localization in HgTe quantum wells: 

E F < m(k F ) 

Let us now consider the limit of low densities, when 
the Fermi energy is smaller than to. This is the limit of 
the conventional Schrodinger equation in the relativistic 
Dirac description. In this case, the spectrum can be ex- 
panded in A 2 k 2 (for definiteness, we consider the upper 
branch) : 



E+(k) = e(k) + y/A 2 k 2 + m 2 (k) 

~ \M\ + Bk 2 + e(k) + ^-. 



(45) 



In this case the pseudospin is almost fixed along the z- 
direction by the "Zeeman field" M, while the off-diagonal 
terms in the Hamiltonian can be viewed as a weak "rel- 
ativistic" correction (analogous to spin-orbit coupling), 
see Appendix |Bj On short scales the system is there- 
fore close to the orthogonal symmetry class and should 
exhibit negative magnetoresistance in sufficiently strong 
magnetic field. 

As discussed above, a similar situation is possible 
for the BHZ Hamiltonian at very high densities, when 
Bkp ^> { | Af | , Akp}, where the spectrum is given by 



E + {k) ~ Bk 2 + e(k) + M 



B ' 



(46) 



However, in realistic systems, when the first term in the 
small- A: expansion becomes dominant, the expansion is 
expected to break down: all higher-order terms are ex- 
pected to be as relevant as the lowest-order one {Bk 2 
here). The analysis of the system at such high ener- 
gies requires a more detailed knowledge of the spectrum. 
Therefore, here we will focus on the controllable case of 
small energies. 

Now we consider the massless Dirac Hamiltonian Ha = 
A(k 

x®x ky<Ty)s z as a perturbation to the diagonal part 
of the BHZ Hamiltonian 



H 



M 



-m(k)a z s z , 



(47) 



dominated by the mass term m(k) = M + Bk 2 . The 
Hamiltonian (|47j) is invariant under the following eight 
TR symmetries: 



Too '■ 


i ^ a Q soO T aos , 


rp'2 

1 oo — 


1, 


(48) 


T 0z : 


O 1 ^ OoSzO 1 (ToS z , 


rp2 

1 0z ~ 


f, 


(49) 


T ■ 

J XX ■ 


O H> cr x Sx0 7 <7 x Sx, 


rp2 

XX 


1, 


(50) 


T ■ 


O H> <7 x s y 7 a x s y , 


rp2 

1 xy — 


-1, 


(51) 


Ty x '■ 


O H-> <7ySxO T (JySx, 


rp'2 

1 yx — 


-1, 


(52) 


Tyy '■ 


O 1— > (JySyO T (JySy, 


rp2 

yy ~ 


1, 


(53) 


T z o ■ 


O i ^ a z soO T a z so, 


rp2 

1 z0 — 


1, 


(54) 


T ■ 

-L ZZ ■ 


O i ^ a z s z O T a z s Zl 


rp2 

ZZ 


1. 


(55) 



Too r 


-L ZZ 1 


Toz r 


° T z o, 


T r 
- 1 - XX 


± VVi 


T r 
J-xy 


•j T 



We see that there are overall six symmetries of the 
"orthogonal" type (yielding WL) and two "symplectic" 
(yielding WAL). However, we have to take into account 
that the Hamiltonian Hm breaks into blocks correspond- 
ing to positive and negative energies (conduction band 
and valence band). The two blocks correspond to eigen- 
values of a z s z equal to ±1. The transport is determined 
only by one of these blocks (where the Fermi energy is 
located). 

Therefore, eight symmetries Eqs. (j48|) - (f55|) split in 
four pairs, 



(56) 



and each pair should be counted only once for the calcu- 
lation of the WL correction since two symmetries of the 
pair become identical when reduced to any of two a z s z 
eigenblocks. This yields 

fo -< 3 - 1 >^'°(T) = -S k (T)' < 57 » 

that is the contribution of an antilocalizing singlet and 
a localizing triplet. This is a conductivity correction for 
two independent copies of orthogonal symmetry class. 

Let us now include the off-diagonal Dirac term Ha = 
A(k x a x + k y cTy)sz. The Hamiltonian Hm + Ha still sepa- 
rates into two blocks with approximate orthogonal sym- 
metry but positive- and negative-energy sectors are now 
mixed. Effective Hamiltonian valid near the Fermi energy 
is derived in Appendix |B] Breaking of the time-reversal 
symmetry within a single block of the BHZ Hamiltonian 
can be traced back to the appearance of an effective mag- 
netic field [with the vector potential given by Eq. (|B7[) ] 
due to the interplay of "relativistic corrections" and dis- 
order scattering. The weak localization correction (1571) is 
modified by the presence of this effective magnetic field, 



T 



(58) 



where 1/ta is the symmetry-breaking rate due to Ha 
(it is analogous to 1/tm introduced above). The value 
of 1/ta can be estimated as follows. Adding disorder 
potential to Hehz, we project the full Hamiltonian on a 
single chiral branch. The off-diagonal part Ha would lead 
to the appearance of the "spin-orbit" impurity scattering 
characterized by 



A 



A ~ - 



I A 2 k 2 F 



t m 2 {kp) 



This estimate yields 



Ak f 



to(/cf) 



(59) 



(60) 
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Alternatively, using the reduced spinless Hamiltonian 
(|B8I) . one can evaluate the averaged phase difference 5 a 
("dephasing action") accumulated due to the effective 
vector potential, Eq. (IB7|) . on the time-reversed paths 
within the path-integral formalism (for dchnitcness, we 
consider the case when the Fermi level is located near the 
bottom of the spectrum, Ef <C \M\): 



/dl ia ( ri ) /dl 2 a(r 2 )\ 

J J I dis 

(J dh x VF(ri) y <fl a x W(r 2 ) 



Using the short-range correlated potential, 

5(rj - r 2 ) 



(61) 



(y(r 1 )y(r 2 )) di! 



one can estimate 



VFk F 
M 2 VT 



2hvt 



(62) 



where ^ ~ |M|/A 2 is the density of states, vf = 
UfA 2 /\M\ is the Fermi velocity, and all the gradients, as 
well as the inverse size of impurity, are replaced by fcp-. 
From S a ~ 1 we find t ~ M 2 vt jvpkp, which yields Eq. 
(fBTij) . Again, the derivation 1 ^ of a Cooperon within the 
kinetic equation approach and diagrammatics confirms 
this estimate. 

The scaling of the corresponding symmetry-breaking 
length Ia = (Dta) 1 ^ 2 with the carrier concentration n 
for energy-independent r has the form: 



I A oc 



(63) 



It is interesting to note that there is a formal relation be- 
tween the symmetry-breaking lengths in the symplectic 
[Eq. (|iD])] and orthogonal [Eq. ([6"3"]l] cases: 



I 



(64) 



In each of the cases only one length is relevant, whereas 
the other is then shorter than the mean free path, imply- 
ing no diffusive dynamics in the corresponding channel. 
Both lengths become of the order of the mean free path 
for Ef ~ m(fcir), where without valley mixing the log- 
arithmic first-order interference correction does not de- 
velop at all (2U phase). 

Let us now include BHZ block-mixing rates t^ 1 and 
TgQ. As we have pointed out above, the eight symme- 
tries ((4*5)) - ([551 of the Hamiltonian Hm split into four 
pairs (|56p . The symmetries within each pair are indistin- 
guishable near the Fermi energy. The general rule for the 
evaluation of the WL correction is as follows. For each 
of the pairs (JSjoJ) , one of the following three situations is 



possible: (i) Both symmetries from the pair are preserved 
(which means that the system splits into two eigenblocks 
of u z s z ). Then the pair contribute to the WL correc- 
tion as a single time-reversal symmetry, (ii) One of two 
symmetries is broken, the other one is preserved. Then 
the situation becomes conventional, and the remaining 
mode gives a usual logarithmic contribution, (iii) Both 
symmetries from the pair are broken. Clearly, there is no 
contribution from this pair to the WL (WAL) correction 
in this case. To summarize, the pair gives a contribu- 
tion as a single Cooperon mode, unless both symmetries 
are broken. The fact that one should break both symme- 
tries to suppress the contribution implies that one should 
add symmetry-breaking times (rather than rates) corre- 
sponding to both symmetries of a pair and then invert 
the result in order to get the mass of the corresponding 
mode. 

In our case additional terms have the following symme- 
tries. The Dirac kinetic term Ha preserves only T xx and 
T xy out of the list pg jHjggj) . The A-term H A = ^o^zS x 
conserves Too, T z q, T xy , and T yy . The Rashba and 
Dresselhaus terms (linear in momentum) conserve T xy 
and T yx . This yields the following result for the WL 
correction: 24 



5a 



e 

2irh 



2 In 



In [ — 

T 4> 



T 

+ — 

TA 

T 
T~SO 



T 
TSO 



TA + T A 



-Jul- 

T0 



(65) 



The only massless mode here again corresponds to the 
T xy symmetry (physical symplectic TR symmetry). Thus 
in the presence of all block-mixing terms, the system be- 
comes a single copy of the symplectic class. Depending 
on the relation between t7 , t^ 1 , and t^q, three pat- 
terns of symmetry breaking are possible: 



(i) 20 — > 2U — > ISp. This scenario is realized when 
t C 7,4 < min[rA, t$o ]■ In this case, at high 
temperatures the system behaves as two copies of 
orthogonal class (20) with doubled WL. In the 
intermediate range of temperatures ta "C t^ <C 
min [ta, T50], the system behaves as two copies of 
unitary class (2U) with no T-dependent interfer- 
ence correction. Finally, at 3> min[rA, tso] the 
two spin blocks are completely mixed and the sys- 
tem reaches its generic spin-orbit state ISp (single 
copy of symplectic class with ordinary WAL). 

(ii) 20 — > ISp. This scenario is realized when 
min[T A ,T5 ] < t a and r < mm[TA,Tso], i-e., 
when the block mixing is faster than the breaking of 
the effective orthogonal TR symmetry within each 
block. As a result, there is no room for the unitary 
class in this case. 

(iii) ISp. When the complete mixing of spin blocks is 
very fast, which happens when tso ^ r : the system 
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FIG. 1: Schematic diagram showing the symmetry patterns in the L — Ef plane. Here L is the system size, phase-breaking 
length L,),, or magnetic length Lh, whichever is shorter. For simplicity, the transport scattering length l tr (vertical dotted line) 
is assumed to be energy independent (this assumption does not affect the "topology" of the diagram). The "phase boundaries" 
are shown by solid curves. Boundaries of the 2U-regions are determined by l m = {Dr m ) 1/2 for 2Sp -> 2U crossover, and by 
I a = (Dt a ) 1/2 for 20 -)■ 2U crossover. Left panel: Inverted band structure (thick quantum well) with no block mixing. Dashed 
line shows energy for which m{kF) = —\M\ + Bkp — 0. Middle panel: Normal band structure (thin quantum well) with no 
block mixing. Right panel: Inverted band structure with block mixing, characterized by Za = (Dta) 1 ^ 2 and lso = {Dtso) 1 ^ 2 
(dash-dotted), assuming ta < tso- The energies Ef = Ei . . . E& (from bottom to top) mark different horizontal cross-sections 
of the "phase-diagram" corresponding to the patterns 20 ISp, 20 — > 2U — ¥ ISp, 2U — ¥ ISp, 2Sp — ¥ 2U — ¥ ISp, and 2Sp 
— ¥ ISp, respectively, that appear with increasing L. The magnetoresistance curves corresponding to these energies are shown 
schematically in Fig. [5] The diagram for the case of normal band structure with the block mixing is qualitatively similar. 




B B B B B 



FIG. 2: Schematic view of the magnetoresistance curves p xx (B) in linear-logarithmic scale. Panels correspond to energies 
Ei . . . E$ from the right panel of Fig. [T] (inverted band structure). The numbers indicate the prefactor a = [nh/ e 2 )da / 8th 
(here th = L 2 H /D) of the logarithmic magnetoresistance in different domains of magnetic field. For the purpose of visualization, 
the crossover regions (where th is of the order of the corresponding symmetry-breaking time) are shown as cusps. 



behaves as a single copy of symplectic class (ISp) 
in the whole diffusive regime (ordinary WAL). 

In the cases (i) and (ii), the magnetoresistance changes 
from positive to negative with increasing magnetic field. 
In the third case, the magnetoresistance is always positive 
in the diffusive regime. 



IV. DISCUSSION AND CONCLUSIONS 

To summarize, we have explored symmetries of the 
Bernevig-Hughes-Zhang Hamiltonian of 2D charge car- 
riers in HgTe/HgCdTe quantum wells in the presence 
of physically relevant symmetry breaking perturbations. 
We have identified regimes of different symmetry in the 
parameter space and evaluated the corresponding quan- 
tum interference corrections to conductivity. Possible 



regimes include 20, 2U, 2Sp, and ISp, with the tem- 
perature dependence of the 2D weak (anti-)localization 
correction given by 5a = a(e 2 /irh) In t>(T), where a = 
— 1,0,1, and 1/2, respectively. These regimes are sum- 
marized in Fig. [T] 

Experimentally, the quantum interference — weak lo- 
calization or antilocalization — shows up most directly in 
the transverse-field magnetoresistance. In particularly 
interesting situations, a symmetry breaking pattern then 
determines a succession of regions of magnetic field with 
different signs and/or prefactors of the weak localization 
magnetoresistance (see Fig. [2]). 

Let us now discuss published experimental dat a 19 i 20 
in context of our findings. In these papers the weak lo- 
calization was studied in structures with both normal 20 
and inverte d 19 ' 20 band gaps. In both cases the authors 
observed weak positive magnetoresistance in low mag- 
netic fields B which could be clearly attributed to WAL. 
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In higher fields, a crossover to negative magnetoresis- 
tance was observed that could be presumably attributed 
to weak localization. This corresponds to the 2 O — > 1 Sp 
symmetry pattern [characteristic for systems with rela- 
tively small carrier concentration, Ep <C m{kp)\ in our 
terminology. The coefficient in front of (e 2 /irh) \tith 
(where tjj — L 2 H /D) was found to be consistent with 
1/2, as expected for the ISp regime. 

More experimental work is clearly needed to explore 
systematically different parameter regimes with different 
types of quantum interference behavior. We hope that 
our paper will stimulate such experimental activity and 
will be helpful in identification of different regimes and 
analysis of quantum interference contributions to conduc- 
tivity. 

When this manuscript was in preparation, we learned 
about preprints Ref. [25ll26l with partly overlapping con- 
tent. 
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Appendix A: Boundary scattering 

In this Appendix we consider the properties of the BHZ 
Hamiltonian with the boundary condition ip = 0, see 
Sec. Ill CI We demonstrate that the symplectic (block- 
wise) time-reversal symmetry is strongly violated by this 
boundary condition. 

Consider the upper block of the BHZ Hamiltonian (to 
minimize the notation we assume A = 1) 



h = 



M - 

k'_r 



Bk 2 



h x thy 
-M + Bk 2 



(Al) 



At energy E = |k| = y/M/B, diagonal elements of the 
Hamiltonian vanish leading to the emergence of symplec- 
tic time-reversal symmetry. There are two plain wave so- 
lutions at this energy with the fixed momentum k x along 
x direction: 



t / J± = ^ 1 J e ik x x±i^/E 2 -kiy^ 



a± = 



k x ± i^/E^ - kl 



E 



(A2) 
(A3) 



These two states are related by the time-reversal opera- 
tion: 



a y ijj±{k x ) = -ia T ip T (~k x ). 



(A4) 



Assume a semi-infinite plane y > with the hard wall 
boundary at y = 0. The eigenstate of the boundary prob- 
lem will contain both the incident wave and reflected 
wave ip+ with some amplitudes. In order to fulfill the 
boundary condition i)\ y =o = 0, we have to add a third 
eigenfunction of the 2D problem exponentially decaying 
in the bulk y > 0. This solution has the form 



1 

Bk 



a 



Jk x x-y/l-B2(E2-kl)y/B 



^1 - B 2 (E 2 - k 2 x ) 



(A5) 
(A6) 



1 + BE 

The scattering state can now be directly constructed: 

(A7) 
(A8) 



ao - a+ 



ao - a+ 



Far from the boundary, the solution ipQ decays and the 
factor r yields the reflection amplitude. 

If the boundary condition preserve time-reversal sym- 
metry, the following state would be another legitimate 
solution of the boundary problem at y — > +oo: 



Applying relations (|A4|) . we obtain 



(A9) 



a y ip* = a y ipZ + r*<T y ip+ 

= -ia + ip + (-k x ) - ia-r*ip-(—k x ) 
a 



ip+(-k x ) 



(A10) 



Now we compare the pref actor in front of k x ) with 
the reflection amplitude at — k x . This will indicate the 
degree of time-reversal symmetry breaking during one 
scattering off the boundary. Actually, the following iden- 
tity holds: 



a+ 



a-r*(k x ) 



-r(-k x ). 



(All) 



The sign here is just opposite to what one would expect 
from the time-reversal invariance. Thus the symmetry is 
maximally violated. 



Appendix B: Relativistic corrections at the bottom 
of the band 



At energies close to the bottom of the band, such that 
E = M + e with e < M, the BHZ Hamiltonian ac- 
quires an approximate orthogonal time-reversal symme- 
try, see Sec. MI Bl It is related to the fact that the spin 
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of the electron is almost aligned with the z axis and very 
weakly depends on the momentum. The Hamiltonian 
can be approximated as h « k 2 /2Af (in this Appendix 
we assume A = 1 and B = 0) up to small relativistic 
corrections. The situation is completely analogous to the 
non-relativistic Shrodinger approximation to the Dirac 
Hamiltonian. Below we derive these relativistic correc- 
tions for the upper block of Eq. @. 

Let us write explicitly the two-component Dirac equa- 
tion with an external potential V: 

(M + V)u + k-V = (M + e)u, (Bl) 
k+u - (M - V)v = (M + e)v. (B2) 

Here k± = k x ±ik y . We now express v from the second 
equation, 



v = (2M + e - V^k+u, 



(B3) 



and use it to obtain a single second order equation for 
the component u, 



[V + fc_(2M + Bk 2 + e - Vy 1 ^] \ 



(B4) 



In order to recast it in the form of an equivalent one- 
component Schrodinger equation, we have to correct the 
normalization of the wave function. The initial spinor 



function is normalized according to \u\ 
rescale the component u as 



u = 1 - 



8AP 



With the help of Eq. (|B3)l, this yields \<j)\ 2 = 1 
We now substitute u from Eq. (|B5[) into Eq. 



1. We 



(B5) 

f 0(fc 4 ). 
Bill and 



expand in powers of k. In order to cancel the terms 
~ ek 2 we have to multiply both sides of Eq. (|B4[) by (1 — 
k 2 /8M 2 ). This finally yields the Schrodinger equation 
hah = e<j) with the Hamiltonian 



h 



2M 



8Af 3 



V 



k-Vk^ 
AM 2 



k 2 V + Vk 2 
8M 2 



(B6) 



First two terms represent the expansion of the relativis- 
tic dispersion relation e = v M 2 + k 2 — M in powers 
of k. The last three terms appear due to external po- 
tential V. In an analogous calculation for the full 4D 
Dirac Hamiltonian, the last two terms of the Hamilto- 
nian would correspond to spin-orbit scattering. In our 
2D case, non-relativistic wave function is just the scalar 
(f> without any spin structure. The "spin-orbit" scatter- 
ing in this case can be rewritten as the interaction with 
a fictitious magnetic field. Let us introduce the vector 
potential according to 



4eM' y 4eM' 
This allows us to represent the Hamiltonian as 

(k - ea) 2 k 4 



(B7) 



h = 

2M 8M 3 
The effective magnetic field 



V- 



2M 



- A. (B8) 
2M 



V 2 F 

_ 4eM 

leads to breaking of the approximate "orthogonal" time- 
reversal symmetry with the rate 1/ta as calculated in 
Sec. IIHB1 see Eqs. (|60|)-([62]>. 
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